Fractal Droplets in Two Dimensional Spin Glasses 
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The two-dimensional Edwards- Anderson model with Gaussian bond distribution is investigated at 
r = with a numerical method. Droplet excitations are directly observed. It turns out that the 
averaged volume of droplets is proportional to l^ with D = 1.80(2) where I is the spanning length of 
droplets, revealing their fractal nature. The exponent characterizing the / dependence of the droplet 
excitation energy is estimated as O'j^ — —0.42(4), clearly different from the stiffness exponent for 
domain wall excitations. 
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^'Droplef is a well-known and useful concept in the 
study of critical phenomena W. The droplet argument 
for spin glasses |0,0] is based on this concept and is one 
of important working hypotheses in the field. Its va- 
lidity is being actively discussed, in particular, for the 
low-temperature phase in three dimensions B . Droplets 
are defined as collective excitations from some pure state 
below the transition temperature. In the case of a typical 
spinglass model in two dimension, however, the system 
is believed to be critical right at the zero-temperature. 
Therefore, it is not clear, a priori, if we can treat droplets 
defined at zero-temperature in two dimensions in the 
same fashion as we treat the droplets below the criti- 
cal point in three dimensions. In two dimensions, most 
of the critical indices can be derived from a single ex- 
ponent Od that characterizes the size-dependence of the 
droplet excitation energy, according to the droplet argu- 
ment or other similar arguments such as the domain wall 
renormalization group theory m. For example, —6d was 
expected to be equal to yt = 1/v while —9s was identified 
with yt pi , where 9s is the stiffness exponent which char- 
acterizes the system size dependence of the domain wall 
excitation energy. A deviation from this standard story 
for two-dimensional spin glass models was first reported 
based on numerical calculation of magnetization at zero- 
temperature with finite field [pi , in which we found that 
the thermal exponent yt does not agree with —9s- Several 
subsequent numerical works jTHlOl] confirmed this result. 
Thus the equivalence among three exponents, yt, —9s and 
—9d, should now be a subject to a close re-examination. 
As for the equivalence between yt and —9d, so far no evi- 
dence of its violation has been found. However, this may 
be simply because identifying droplets is technically so 
difficult that nobody has actually succeeded in direct ob- 
servation of droplets large enough to discuss asymptotic 
behavior. In fact, if one adopts the original definition g] 
of droplets and apply a simple combinatorial optimiza- 
tion algorithm, it would be impossible to obtain droplets 
larger than a few lattice spacings. 

In this letter, we focus our attention on the EA model 



on a square lattice 
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where Jy 's are quenched Gaussian random variables with 
the mean value of zero and the standard deviation of J . 
We take J as the unit of energy throughout the present 
letter. It should be noted that the one-parameter scal- 
ing which was derived from various "pictures" such as 
the droplet picture and the domain wall renormalization 
group argument can be also derived within the framework 
of standard finite size scaling theory without any assump- 
tions or pictures. We start with the following standard 
form for the partition function. 
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where Z, T, H, L are the singular part of the partition 
function, temperature, magnetic field, and system size, 
respectively. We hold (0) to be a defining equation of 
the exponent yt- One of the two scaling exponents can 
be fixed by using the fact that the magnetization at zero 
temperature and zero magnetic field is proportional to 
ljd/2 because of the absence of non-trivial degeneracy. 
This leads to yh = yt + d/2. In previous works, we mea- 
sured the magnetization M{T = 0, H, L) and the spin- 
glass susceptibility xsg{T,H = 0, L) [Q which are both 
derivatives of log Z with respect to T and H multiplied by 
some known factors. By matching the numerical results 
with the asymptotic form derived from (yj), we obtained 
yt ~ 0.48 in both calculations. 

On the other hand, it is often assumed, or for some 
models proved, that the excitation energy Ew of domain 
walls induced by twisting the system has the same scaling 
behavior as the singular part of the free energy. This 
means 9s = —yt where ^5 is defined by 

Ew{L)ccL^'' {T = 0). 

In the early stage of the study on the present issue, the 
domain wall energy was measured for various sizes p] 
with the hope of obtaining an estimate of yt ■ The stiffness 



exponent was estimated as —0s — 0.29 M. This result 
was reconfirmed later by other groups |^,]ll| with larger 
scale computations. Therefore, it is now rather clear that 
—0s and yt are two different exponents. 

Fisher and Huse defined ^ a droplet of the scale A 
including a given site ? as a cluster of spins (with i among 
them) with the smallest excitation energy that contains 
more than A'' and less than (2A)'' spins. The basis of the 
argument then is the following scaling form that describes 
the excitation-energy distribution of droplets of the scale 
A: 



Px{Ex) = 



TA»" V TA«^ 



(2) 



where T is some constant, 0o is the droplet exponent, 
and P{x) is the scaling function which is continuous and 
non-vanishing at a; = 0. Since the droplet "size", A, is 
defined to be proportional to (volume)"'^/'', it is not neces- 
sarily a spanning length of droplets because the volume of 
a droplet may in general have a non-trivial fractal dimen- 
sion. In fact, in what follows we demonstrate that this 
is the case, i.e., a typical large droplet in two dimensions 
occupies only an infinitesimal fraction of the volume of 
the box that contains it. Therefore, for a given droplet, 
we have at least two different length scales, A and the 
spanning length. We refer hereafter to the fractal di- 
mension of the droplet volume as D. Then, we have the 
relation 



V = X'^ozl° 



where I is the spanning length of the droplet. 

Since the above-mentioned definition of droplets is in- 
convenient from computational point of view, we have 
adopted the following alternative definition. First we 
consider an L x i system with free boundary condition. 
We define a droplet of scale L as the cluster of spins that 
has the smallest excitation energy among those which 
contain the central spin and contains no spins on the 
boundary. 

If we assume the conventional droplet argument, the 
central spin is surrounded by droplets (in the original 
definition) of various scales. Since perimeters of smaller 
droplets are far from the boundary of the system they 
must be relatively free from the influence of the con- 
straint imposed on the boundary. Therefore, smaller 
droplets in a finite system should have the same scal- 
ing properties as they would in an infinite system. On 
the other hand, the fixed boundary imposes strong re- 
strictions upon the shapes that larger droplets can take. 
As a result, they have larger excitation energy than they 
would in an infinite non-restricted system. Therefore, a 
crossover size of droplets must exist and it separates the 
region of smaller sizes in which one observes a correct 
asymptotic behavior from the other region of larger sizes 
that is strongly affected by the boundary. 
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FIG. 1. 16 randomly chosen droplets. 



In the spirit of effective mutual independence among 
the droplets of different scales discussed in H , it is nat- 
ural to expect that the crossover size corresponds to the 
largest droplet that can fit into a box of L x L, i.e., 
the droplet whose spanning length is of the order of L. 
Then, considering the fact that larger droplets tend to 
have smaller excitation energy, one may also expect that 
the spanning length of the droplet identified by the pro- 
cedure described above would be proportional to L. We 
here define 6*^ by the following L dependence of the av- 
erage excitation energy of the droplets in our definition, 

[Ed{L)] oc L^'" , 

where [• • •] stands for the average over samples. Compar- 
ing this to the original definitions of the droplet size and 
droplet exponent, we obtain 



W 



(3) 



because L is proportional to the spanning length of the 
droplet. 

In order to observe droplets in the new definition, we 
first compute the ground state with free boundary con- 
dition, and take it as the reference spin configuration. 
Then we compute the ground state with the constraint 
that the spins on the boundary are to be fixed as they 
are in the reference state while the central spin is to be 
fixed opposite, thereby forcing a cluster of spins includ- 
ing the central spin to fiip. For a system with the free 
boundary condition, polynomial-time optimization algo- 
rithms are available whereas for the systems with con- 
straints no such algorithm is known. In fact, the two 
dimensional spin glass problem with general constraints 



has been proven to be NP hard |1^ . Therefore, we have 
employed the rephea optimization M , which is a heuristic 
optimization algorithm based on the idea of renormaliza- 
tion group. The details of the algorithm are described 
elsewhere P,p^. 

Since the algorithm is only heuristic, we need to esti- 
mate the rate of obtaining false ground states and how 
much such states may affect the average values of vari- 
ous quantities of interest. For this purpose, we compared 
our results with exact solutions when they are available. 
We have actually found no error with the settings of al- 
gorithmic parameters that are used in the computation 
presented in this letter. We detected some errors only 
when we loosened the termination condition. Even in 
that case, the rate and the amount of the errors turned 
out to be so small that their effect is well within the final 
statistical errors. Therefore, we believe that the errors 
are negligible even for the largest systems (L = 49), for 
which no exact result is available for comparison. 

In Fig.y, we present 16 randomly chosen droplets in the 
system of size L = 49. We can see that the size of the 
droplets varies depending on the sample. Some of them 
contains only one spin. This reflects the fact that |6'd| 
is small and the excitation energy does not very strongly 
depend on the droplet size. Frequency of having small 
droplets decreases as a function of the system size. On 
the other hand, we also see that larger droplets whose size 
is comparable even to the system itself occupy a consid- 
erable fraction in the whole population of droplets. In 
these larger droplets, many fingers and overhangs can be 
observed, which already suggests the fractal nature of the 
droplets. 

In order to check that the spanning length of a typical 
droplet is really proportional to the system size as we 
expect, we compute two length scales Iq and h- The 
length Iq is the spanning length itself, i.e., the difference 
in x-coordinates of the left-most site and that of the right- 
most site in a droplet, assuming that the a;-axis stretches 
from left to right. On the other hand, the length I2 is 
defined as 



h^^{ix-{x)y}. 



Here, 
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where Xi is the x-coordinate of the site i. In FigJ^, these 
two length scales averaged over many samples are shown 
as a function of the system size L. A correction to scaling 
can be observed for small system sizes L < 17. In partic- 
ular the slope for ^0 niay seem slightly larger than unity 
at a first glance. However, since the slope cannot asymp- 
totically exceed unity for an obvious reason, we consider 
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FIG. 2. Various length scales of droplets. The statistical 
errors for L = 49 are comparable to the symbol size whereas 
those for other data points are much smaller than the symbols. 



that /o is proportional to L in the large L limit. The lin- 
earity for I2 is much better and indicates proportionality 
to L. 

The averaged length, P, of the boundary of the droplet 
and the averaged volume, V are also plotted in Fig.g. The 
linearlity in the logarithmic scale is rather good for both 
P and T/^/^ and both of them have slopes different from 
unity. For P, we obtain P oc L^' with the surface fractal 
dimension 

D, = 1.10(2). 

For 1/, we estimate the fractal dimension as 

D = 1.80(2). (4) 

The fractal dimension D of droplets is certainly smaller 
than d — 2. Thus, we conclude that the droplets at the 
critical point T = have fractal nature in the volume as 
well as in the perimeter. 

Finally, we measure the droplet excitation energy 
Ed(V). Similar to the size and the shape, the droplet ex- 
citation energy has a broad distribution. When rescaled 
with the average value Ed (L) , the histograms of excita- 
tion energies for various system sizes fit on top of each 
other, showing the validity of the form (H) (See Fig.^). In 
addition, we observe that the scaling function P{X) has 
a non-vanishing value at X = 0, satisfying a necessary 
condition for the droplet argument to be valid. 

Since the average droplet excitation energy does not 
strongly depend on the system size, it is subject to rel- 
atively large correction to scaling and the estimate of 
9'jj can be affected by the deviation from the asymp- 
totics in small L region. Therefore, the estimate depends 
on which we choose among the above-mentioned length 
scales that differ from each other for small sizes. When 




4.0 



1.0 2.0 3.0 

X = E,(L)/[E,(L)] 

FIG. 3. The distribution of droplet excitation energies 
rescaled with the average values. 
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FIG. 4. The average droplet excitation energy. 



log Ed is plotted against log L the slope is estimated to be 
d'jj ^ —0.42, whereas it becomes —0.38 and —0.45 when 
^0 and I2, respectively, are chosen instead of L. These 
three plottings are shown in Fig.Q together with straight 
lines with slopes of corresponding estimates of 6'jy . From 
these results we quote the following value as the estimate 
of the droplet exponent 9'jj: 

- ^i, - 0.42(4). (5) 

Assuming (m and using the value of D in (0) , we obtain 
the following estimate of 9jj , 

- 0i3 = 0.47(5) (6) 

in good agreement with previous estimates of yt such as 
yt = 0.48(1)0. 

To summarize, we have generated droplets in two di- 
mensions for the EA model with Gaussian bond distri- 



bution. We have found the droplets non-compact in con- 
trast to what is assumed in the conventional droplet ar- 
gument for higher dimensions. The exponent —O'j-, has 
been estimated to be slightly smaller than the thermal 
exponent in the absolute value, and this difference can 
be well understood within the framework of the finite- 
size scaling and the standard droplet argument, by taking 
into account the fractal nature of droplets. It is rather 
surprising that the droplet argument seems to be valid 
for two dimensional critical behavior in spite of the qual- 
itative difference between fractal droplets at the criti- 
cal point shown in the present letter and the off-critical, 
presumably compact, droplets in higher dimensions. The 
disagreement between the droplet exponent and the stiff- 
ness exponent also remains puzzling. At least, however, 
the present result suggests that a domain wall oi a L x L 
system may not necessarily be considered as an object 
of scale L because of the existence of different ways for 
measuring the size of the same object. In other words, 
at present we do not know which way of measuring we 
should use for comparing the size of a domain wall and 
that of a droplet. 
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